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ABSTRACT 

The vacuum-polarization corrections to rates for nuclear-fusion reactions in 
the pp chain and in the CNO cycle are calculated. For the reactions of particu- 
lar importance to the solar-neutrino problem, the "^He(^He, 2p)'^He, "^He(a, 7)^Be, 
'''Be(p, 7)^B, and ^^N(p, 7)^^0 reactions, we find the magnitude of the effect to be 
less than 2%. The effect of VP on all the other reaction rates is expected to be 
of a similar order of magnitude. We discuss how these results affect the predicted 
fluxes of solar neutrinos. 
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The cross sections for the nuclear reactions in the pp chain and in the CNO 
cycle are important for stellar-evolution calculations of main-sequence stars [1] 
and, in particular, for the calculated flux of solar neutrinos [2]. The energies at 
which these cross sections are measured in the laboratory are generally higher 
than the energies of interest in the Sun. The cross sections at solar energies are 
determined by extrapolating the measurements to lower energies using the Gamow 
penetration factor, which expresses the probability of quantum tunneling through 
the Coulomb barrier. In this paper, we calculate the effect of vacuum-polarization 
(VP) corrections to the electrostatic potential on the rates for nuclear reactions in 
the pp chain and in the CNO cycle. We show that such corrections lead to only 
small changes in the predicted fiux of solar neutrinos. 

In addition to the Coulomb potential, Vq = e^/r, there is an additional contri- 
bution to the electrostatic potential, the Uehling potential, Vypir) [3], which arises 
from quantum corrections. The complete electrostatic potential for two nuclei with 
charges Zi and Z2 separated by a distance r is 
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where a is the fine-structure constant, and 
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where mg is the electron mass. The function 7(r) has the limiting forms 

/(r) = — 7 — 5/6 — In(mer), for mgr <^ 1, (3) 

and 

The function 7(r) has a logarithmic singularity for very small radii, and exhibits an 
exponential falloff (arising from the exchange of a virtual electron-positron pair) 
for r > l/2me. The probability for tunneling through the electrostatic barrier is 
affected by the presence of the Uehling potential [4] . 
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The energy dependence of a non-resonant fusion cross section is ordinarily 
written (see, e.g., Refs. [2] and [5]), 



aiE) ^ ^ri^\E), (5) 



where S{E) — S{0) + S'{0)E is a slowly varying function of E, and S{0) and 
S'{0) are determined by fits to experimental measurements of the cross section. 
The quantity r^^\E) is the usual Gamow penetration factor, the probabihty of 
tunneling through the Coulomb barrier. In practice, S{E) is deduced by laboratory 
measurements of the cross section at energies of order 100 keV to several MeV, 
and the cross section is then extrapolated to energies, 0(10 keV), typical of solar 
reactions, through Eq. (5). In this paper, we use a WKB approximation similar 
to that used by Gould [4] to estimate the vacuum-polarization corrections to the 
Gamow penetration factor for any binary reaction. 

For a pure Coulomb potential, the Gamow penetration factor for reacting nuclei 
with center-of-mass energy E is simply 



r(o)(£;) = exp 
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where h — ZxZ^g^ jE is the turning point radius, which is determined by Vc{^ = 
E. Here, rj = ZiZ2e'^ /hv, where v = {2E/fi)^^'^ is the relative velocity of the 
incoming nuclei. Also, /j, — MpAiA2/ {Ai +A2) is the reduced mass of the system, 
and Ai and A2 are the atomic mass numbers of the reacting nuclei. 

To include vacuum polarization in the cross section, we make the substitution 

r(o)(E) ^r(E), (7) 

in Eq. (5), where T{E) is the probability of tunnehng through the complete elec- 
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trostatic potential. Now, 



r{E) = exp 



bvp 

2 

h 



J [2ij{Vc + Vvp-Ep''dr 



r(o)(^)[i-A(^)], (8) 



where the turning-point radius is now given by Vc{bvp) + Vvpibvp) — E, and 
can be written hyp — h[l + 2Q;/(6)/37r]. We expand the integral in powers of the 
fine-structure constant, a, and use the result for hyp to find 



A(£;) = ^ / -^^dx. (9) 





The function A(E) is plotted in Fig. 1 for the ^He("^He, 2j9)^He (solid line), ^He(Q;, 7)^Be 
(short-dash line), "^66(^,7)^6 (long-dash line), and ^^'N{p,'j)^^0 (dot-dash line) 
reactions. The function A(i?) initially increases from its value at ii^ = with in- 
creasing energy, reaches a maximum at an energy Emax and then decreases roughly 
as E~^/'^ at higher energies. Although the lower limit of A(i?) is very conserva- 
tively given by (at very large energies), the lower limit of A(i?) for energies at 
which measurements are performed is actually not much smaller than A{Emax)- 

The most probable energy of interaction for nuclei in the core of the Sun is [2] 

Eq = 1.22 [Zfzi{f,/mp)Ti]^/^ keV, (10) 

where Tq is the temperature in units of 10^ K. We take Tq — 14 (the temperature at 
which energy production is maximized). To self-consistently determine the effect 
of vacuum polarization on reaction rates in the Sun, the VP correction must be 
included in the analysis of the data from which the low-energy cross sections are 
extrapolated. This is done by fitting the data to Eq. (5) using r instead of r(°). The 
cross section is then evaluated at energies typical of nuclear reactions in the Sun 
with the S{E) obtained and with r(£^) rather than r(°)(£'). We then determine 
the effect of VP by comparing this result for the cross section with the standard 
result. 
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For the ^He('^He, 2j9)^He reaction, we used data from Ref. [6] and found that 
including VP self-consistently in the entire analysis decreases the reaction rate by 
about 0.2%. For the ^He(a, 7)^Be reaction, we used the Parker and Kavanagh data 
set [7] and found that VP decreases the reaction rate by 1.6%.* The decrease is due 
to the fact that the data set is weighted at higher energies where the VP correction 
is smaller than that at solar energies. The magnitude of our result is also much 
smaller than 4.6%, which was obtained in Ref. [4]. The discrepancy is due to the 
fact that we included VP in fitting the data, whereas Gould evaluated the correction 
to the cross section at solar energies without correcting for the effect of VP on the 
data. In addition, our value of A{Eq) = 5.3% differs from Gould's value (4.6%) 
because Gould evaluated the correction to the velocity factor in the denominator 
of the WKB approximation to the wave function. In the standard treatment, 
the Gamow penetration factor includes only the leading exponential factor in the 
WKB approximation. Therefore, for consistency with the standard experimental 
procedure, we have not included the velocity factor. For the ^Be(p, 7)^B reaction, 
we used the data set of Kavanagh et al. [8], and found that VP decreased the 
reaction rate by 0.1%. In order to obtain this result, we used only data points at 
energies below the resonance. This is consistent with the treatment in the most 
recent analysis of the ^Be(p, 7)^B S(Q) factor [9]. The ■'^^N(p, 7)^^0 reaction is the 
slowest process in the CNO cycle, so the rate for production of neutrinos from 
^^N and ^^O decays is controlled by the rate for this reaction. The best-estimate 
of the astrophysical S factor for this reaction comes from a data set made up of 
measurements from 0.2 to 3.6 MeV [10]. We used a simulated data set made up 
of 18 uniformly spaced data points spread over this energy range and found that 
VP decreased the solar reaction rate by 0.8%. For each of these four reactions, the 
effect of VP is much smaller than the uncertainty in the measured value of S'(O). 

Our data sets differ slightly from those used to obtain the current best-estimates 
for the low-energy cross-section factors. Our results on the effect of VP would be 
altered by negligible amounts if the exact data sets are used. 

For the other reactions in the pp chain and in the CNO cycle, the effect of 
including VP in the entire analysis can be crudely estimated by comparing the 

★ We are grateful to P. Parker for providing the ^He(ci!,7)''Be and ^Be(p, 7)^B data sets. 
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relative magnitudes of the VP correction, A(£'), at solar energies to that at energies 
at which the measurements are performed. Of course, measurements are performed 
over some range of energies, and the VP correction may change considerably over 
this range, but conservative upper and lower bounds on the VP correction can 
be provided. We do so by comparing A{E) at solar energies with conservative 
estimates of the maximum and minimum VP corrections to the data points. 

In Table 1, we list the VP correction to the Gamow penetration factor at solar 
energies, A{Eq), the most probable energy of interaction, £^0, the maximum of 
A{Emax) = ^{E), and the energy, Emax at which this occurs. In addition, we have 
listed the best-estimates, Syp, of the effect of VP on the rates for the reactions 
that we have analyzed more carefully. In all cases, A{Emax) is only slightly larger 
than A(£'0), and E^ax > Eq. For each reaction, the effect of VP decreases the 
reaction rate by no more than A{Eq) (which is always much less than 6%), and 
may increase the reaction rate by no more than A{Emax) — ^{Eq) (which is always 
less than 1%). The effect of inclusion of VP in the entire analysis should always be 
much smaller than A(£^0). This is because A{E„iax) usually occurs near the lowest 
energies at which measurements are performed, and A{E) ~ E~^/'^ remains non- 
negligible even at higher energies (see Fig. 1). This is illustrated by the results for 
the reactions that we have studied more carefully. Therefore, VP should generally 
have no more than an 0{1%) effect on nuclear-reaction rates in the Sun. 

The rate for the initial — > '^R + e'^ + Ue reaction is also affected by vacuum 
polarization [4,11]. The cross section is too small for this reaction to be observed 
in the laboratory, so the rate for the solar pp reaction is calculated instead of 
extrapolated from measurements at higher energies. Therefore, the effect of VP on 
the pp rate is determined in a different manner, and the result is that VP decreases 
the rate by 0.6% [11]. We hst this best-estimate, as well A(£'0) and A{Emax), in 
the Table. 

Our results imply that VP will have little effect on the predicted flux of solar 
neutrinos. The flux of neutrinos from ®B decay depends most sensitively on the 
nuclear-reaction rates. The dependence of the neutrino flux may be written [2] 

0(«B) oc S^,'-%,'-''SltSh^ (11) 
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where Sn, iSss, 1534, and Sn are the low-energy cross-section factors for the pp, 
^He('^He, 2p)'*He, "^He(a, 7)'^Be, and ^Be(p, 7)^6 reactions, respectively. Inserting 
our best-estimates for the effect of VP on the low-energy cross-section factors, we 
find that inclusion of VP in the analysis increases the predicted flux of ®B neutrinos 
by only 0.2%. 

The dependence of the neutrino fluxes from the other reactions in the pp chain 
and CNO cycle on small changes in the low-energy cross section factors are simi- 
larly determined [2]. The flux of ^Be neutrinos will be decreased by 0.7% by VP 
corrections to the nuclear-reaction rates; the neutrino fluxes from ^^N and ''^^O 
decay will be increased by 0.8% and 1.0%, respectively; and the flux of pp neu- 
trinos will be increased by 0.02%. These result in an increase of 0.01 SNU in the 
predicted event rate for the chlorine experiment and a decrease of 0.11 SNU in the 
predicted event rate for the gallium experiments. 

In summary, we have evaluated the effect of vacuum polarization on the deter- 
mination of nuclear-fusion reaction rates in the Sun. We find that the magnitude of 
the VP correction is never more than 2% for the reactions that are most important 
for solar-neutrino calculations. The VP correction has only a small effect on the 
calculated solar-neutrino fiuxes. 

We acknowledge useful conversations with F. Dyson and R. J. Gould. We are 
grateful to P. Parker for providing several useful data sets. We are also grateful to 
the referee for pointing out an error in an earlier version of the manuscript. MK 
was supported by the Texas National Laboratory Research Commission, and by 
the DOE through Grant No. DE-FG02-90ER40542. JNB was supported by the 
NSF through Grant No. PHY92-45317. 
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FIGURE CAPTIONS 

1) Plot of A(£^) versus the center-of-mass energy. The sohd hne is for the 
^He("^He, 2p)^He reaction, the short-dash line is for the ^He(Q;, 7)^Be reaction, 
the long-dash hne is for the ^Be(]9, 7)^B reaction, and the dot-dash hne is for 
the ^^N(p, 7)1^0 reaction. 
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VP Corrections to Nuclear-Reaction Rates 


Reaction 




A(E0) 


Emax 


^{Emax) 


^VP 


p + p^'^H + e+ + iye 


6 


1.4 


11 


1.5 


-0.6 


^He +^ He ^ a + 2p 


20 


5.0 


46 


5.0 


-0.2 


^He +4 He Be + 7 


21 


5.3 


46 


5.4 


-1.6 


'Li + p^ 2a 


14 


3.3 


31 


3.3 




^Be + p B + 7 


17 


3.8 


46 


3.8 


-0.1 


^He + p He + e+ + i/g 


10 


2.5 


21 


2.5 




'^C(p,7)^^N 


23 


4.8 


66 


4.8 




'^C(p,7)^4p^ 


23 


4.8 


66 


4.9 






25 


5.2 


76 


5.3 


-0.8 


^^N(p,7)^'^0 


25 


5.2 


76 


5.3 






25 


5.2 


76 


5.3 




^^0(p,7)^^0 


28 


5.5 


91 


5.6 





Table 1. Vacuum-polarization corrections to nuclear-reaction rates. £^0 and Emax are 
in keV, and A.{E) and 5vp are in %. 
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